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ABSTRACT 

The number of well-studied early-type pre-main-sequence objects is very limited, 
hampering the study of massive star formation from an observational point of view. 
Here, we present the results of VLT/FORS2 spectropolarimetric and VLT/X-shooter 
spectroscopic observations of two recently recognised candidate Herbig Be stars, PDS 
27 and PDS 37. Through analysis of spectral lines and photometry, we find that these 
two objects are hot, 17500 ± 3500 K, have large radii, 17.0 ± 4.0 and 25.8 ± 5.0 R 0 , 
and are very massive, 15.3 (+5.4, -4.4) and 21.1 (+11.0, -5.3) M 0 for PDS 27 and PDS 
37, respectively. This suggests that these two objects are very young in their evolu¬ 
tion and may become O-type stars. Their youth is supported by their high accretion 
rates of the order of 1CU 3 -10 -4 ' 5 M 0 /yr. A change in linear polarisation across the 
absorption component of Ha is detected in both objects. This change indicates that 
the circumstellar environment close to the star, at scales of several stellar radii, has 
a flattened structure, which we identify as an inner accretion disc. Strong variability 
is seen in both objects in many lines as further indication of an active circumstellar 
environment. 

Key words: techniques: polarimetric - techniques: spectroscopic - circumstellar 
matter - stars: pre-main-sequence - individual: PDS 27 individual: PDS 37. 


1 INTRODUCTION 


Herbig Ae/Be (HAeBe) stars are optically visible pre-main- 
sequence (PMS) stars with masses roughly between 2 and 
10 Mq . With their intermediate masses, they bridge the gap 
between low-mass stars whose formation is well understood 
and high-mass stars who se formation is poorly understood 
(Z innecker fe Yo rkc 20071 1. HAeBe objects were first identi¬ 
fied by Herbid ( 1960T) as stars of spectral type A or B with 
emission lines and which illuminate a bright nebula in their 
surroundings. Their spectral energy distributions (SEDs) are 
characterised by an infrare d excess due to the dust in the 
circumstellar environment dWaters fe Waelkenslll998l ). 

HAeBe objects not only have some commo n character¬ 
istics with high-mass objects such as clustering (|Testi et al.1 
11999l l but also have spectro polarimetric chara cteristics in 
common with T Tauri stars dVink et al.l 1200^ 1. Therefore, 
the study of the formation of HAeBe stars can bring out the 
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differences in the mechanisms of the formation of low- and 
high-mass stars. In addition, because they are closer and less 
embedded than massive young stars, HAeBe stars can act 
as a powerful test bed for the formation of massive stars. 

To study the star formation mechanism as a function 
of mass, it is necessary to sample a large number of HAeBe 
stars, with many representatives in all mass bins. The k nown 
sampl e of HAeBe stars [see e.g. the catalogue bv lThe et al.l 
(1994)] contains many A-type and late B-type objects, but 
there is a dearth of the earlie st type, most mas s ive H erbig 
Be (HBe) stars (also see e.g. lOchsendorf et ahl d201ll )). In 
this paper we present an in-depth investigation of two early 
HBe star candidates identified by the Pico dos Dias Survey 
(PDS). Based on their observational characteristics, PDS 
27 (also known as DW CMa, RA 07:19:36, Dec. —17:39:18) 
and PDS 37 (also known as Hen 3-373, RA 10:10:00, Dec. 
—57: 02:07.3) were pro posed to be early type Herbig stars 
by IVieira et al.l d2QQ3|). The young n ature of PDS 27 was 
also suggested by lSuarez et al.1 (|2006l l. Both stars were also 
classified as young in the Red MSX Source ( RMS) survey for 
Mass ive Young Stellar Obje c ts (M YSOs) (|L umS( len et al.1 
120131 ). In 2011, IVieira et al.1 d201ll ) considered PDS 27 to 
be an evolved star instead, but ruled out an evolved nature 
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for PDS 37. However, as we will demonstrate in this paper, 
both objects are very similar, and we proceed under the 
assumption that they are young stars. 

Using optical photometry, the temperature has been es¬ 
timat e d to be around 220 00 K for both objects (IVieira et al.l 
l2003l f. Ivieira et al.1 d200Hh also estimated their distance and 
it was found to be 1100 pc for PDS 27 and 720 pc for PDS 
37; PDS 27 is associated with CMa st ar-forming region while 
PDS 37 is associated with C 282.4+0. Illee et all ( 2013 1 stud¬ 
ied the near-infrared (NIR) first overtone CO emission of 
PDS 37, and assuming that it originates from the circum- 
stellar disc, they f ound the disc to be hig hly inclined to the 
line of sight (80° 1. [Rodrigues et al.l J2009l l determined an in¬ 
trinsic polarisation of 2.77% for PDS 37, suggesting that the 
circumstcllar environment is not spherical. 

Here we present a spectroscopic and spectropolarimetric 
study of these two objects. We will use the large wavelength 
coverage provided by X-shooter to determine the stellar 
parameters, evolutionary status and accretion rates. Spec- 
tropolarimetry is a very powerful technique that can probe 
very small unresolved spatial scales of the order of stellar 
radii and allows one to extrac t information on th e geom ¬ 
etry of thi s mate r ial (s ee e.g. lOudmaiier fe Drew! dl999ll : 
IVink et al.1 (12001 l2005lf f. A new element in the present 
study is that the medium-resolution spectropolarimetry has 
a wavelength coverage from 4560 to 9480 A, which is much 
larger than presented in any previous study at similar spec¬ 
tral resolution. 

This paper is organised as follows: in Section 2 the ob¬ 
servations and data reduction are discussed. The results are 
presented in Section 3. Section 4 presents an analysis of the 
data followed by a discussion in Section 5. Finally, conclu¬ 
sions are drawn in Section 6. 


2 OBSERVATIONS 
2.1 Spectropolarimetry 

The data of PDS 27 and PDS 37 were taken as part of a 
wider spectropolarimetric investigation whose results will be 
presented in a future paper (Abakakr et al., in preparation). 
The targets were selected based on their strong emission 
lines, allowing an investigation into their spectropolarimetric 
behaviour. 

The spectropolarimetric data for our targets were ob¬ 
tained during the nights of 2012 February 4 and 5 with 
the FORS2 spectrograph mounted on ESO’s VLT in Chile. 
Three gratings were used with a 4096 x 2048 pixel CCD to 
sample the entire optical range. The 1400V grating with cen¬ 
tral wavelength 5200 A covered the spectral range from 4560 
to 5860 A and includes H/3. The 1200R grating was used for 
the spectral range from 5750 to 7310 A, with central wave¬ 
length 6500 A, and includes notable lines: Ha, [O I] 6300 A, 
and He I lines at 5875, 6678 and 7065 A. Finally, the 1027Z 
grating covers 7730 - 9480 A, at central wavelength 8600 A, 
and this includes the calcium triplet lines and the hydrogen 
Paschen series. A 0.5 arcsec slit was used in all cases, which 
resulted in a spectral resolution of R = 3800, 4800, 4800 in 
the V, R and Z bands respectively. 

FORS2 is equipped with an optical polarimeter in or¬ 
der to analyse the linearly polarised components. The po¬ 
larisation optics consist of a rotating half-wave plate and a 


calcite block. The calcite block separates the incoming light 
into two perpendicularly polarised light beams, the ordinary 
(O) and extraordinary (E) beam. The half-wave plate is ro¬ 
tated to measure polarisation at four different angles: 0°, 
45°, 22.5°, and 67.5°. One complete set of observations con¬ 
sists of four exposures, one at each of the four angles. To 
avoid saturation in the bright emission lines, several sets of 
observations were taken with short exposure times. 

The polarisation accuracy is determined only by pho¬ 
ton statistics; in order to reach a relative accuracy of 0.1%, 
one needs to detect roughly one million photons per resolu¬ 
tion element. The targets are relatively faint, with V-band 
magnitudes of 13.0 and 13.5 for PDS 27 and PDS 37 re¬ 
spectively d Vieira et al.ll2003l ). so the total integration times 
were comparatively long with 1.6 and 2.4h total. For PDS 
27 the exposure times are broken down into 8 x 340s for the 
V-grating, 24x50s for the R-grating and 16x120s for the 
Z-grating. Similarly for PDS 37 the times are 16x340s for 
the V-grating, 24 x 70s for the R-grating and 16 x 100s for the 
Z-grating. Polarised and unpolarised standard stars were ob¬ 
served throughout the night. The exposure times for these 
were much shorter as it was only the continuum values that 
were of interest. 

The spectropolarimetric data reduction was carried out 
using gasgano and FORS pipeline v4.9.18. We have used 
two recipes: one dedicated to the calibration and the other to 
extract the science spectrum. The reduction procedure con¬ 
sists of bias subtraction, flat fielding, extraction and wave¬ 
length calibration of the O and E spectra. The Stokes pa¬ 
rameters, obtained from the O and E beams, lead to the 
percentage of linear polarisation, P, and polarisation angle, 
6, according to the following equations: 

P=\fQ 2 + U 2 (1) 

6 = \ arctan (2) 


Finally, the spectropolarimetric data were imported 
into the iraf and polmap package for measurement and 
analysis. 

FORS2 is mounted on the UT1 Cassegrain focus, re¬ 
sulting in a relatively stable instrumental polarisation. To 
improve signal-to-noise and the polarisation accuracy, the 
individual data sets were combined. The instrumental po¬ 
larisation is determined from the observation of unpolarised 
standard stars. This was found to be ~ 0.16%. The polari¬ 
sation angle has an absolute error of ~ 0.5° as determined 
from observing polarised standard stars. We did not per¬ 
form a correction for instrumental polarisation as it is small 
compared to the continuum values of the targets, while our 
main emphasis lies on the differential polarisation between 
lines and continuum. 


2.2 Spectroscopy 

Spectra of PDS 27 and PDS 37 were obtained using the 
medium-resol ution spectrograph X-shooter, mounted on the 
VLT, C hile (iV ernet et abl 201 ill as part of a larger pro¬ 
gramme (lOudmaiier et alll2011 : lFarilamb et al.ll~2015lb One 
of the main strengths of the instrument is its huge wave¬ 
length coverage, taken simultaneously, over three arms: 3000 
- 5900A, UVB; 5300A - 1.0/im, VIS; and 1.0 - 2.4/rm, NIR. 
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Figure 1 . A selection of observed spectral lines from the X-shooter data is shown for both targets (for PDS 27 only the X-shooter(b) 
data set is shown for clarity). The spectra are shifted such that their systemic velocity (see the text) is 0 kms -1 . The red solid line denotes 
PDS 27, while the blue dashed line is P D S 37. The black dash-dotted line is a Kurucz-Castelli (KC) model atmosphere for a 21000 K 
type star llKurucdl 19931 : [Castclli <V Kuruc3l2004 ~): this is an upper limit to the temperature (this limit will be discussed in Section 4.1.3). 
The spectra of both objects are normalised to one but PDS 27 is shifted by an arbitrary value for clarity. Note the similarity of both 
objects’ spectra, with a difference being the number of double-peaked lines of Fe II in PDS 37. Blueshifted absorption reaching -700 and 
-400 kms -1 is observed in Ha and helium for PDS 27 and PDS 37, respectively. 


The smallest slit widths available (0.5 arcsec - UYB, 0.4 arc- 
sec - VIS and 0.4 arcsec - NIR) were used and provided the 
highest spectral resolution available for X-shooter: R=9100, 
17400 and 11300, respectively. Both objects were observed 
in nodding mode in an ABBA sequence. Two observations 
were made of PDS 27, on the evenings of 2009 December 
18 and 2010 February 24 (these data sets will be referred 
to as X-shooter(a) and X-shooter(b), respectively, for the 
remainder of this paper). The exposures were the same for 
each observation with a breakdown across the arms of UVB 


- 300sx4, VIS - 300sx4 and NIR. - (2sx6)x4. PDS 37 was 
observed on the evening of 2010 March 31 with exposures of 
UVB - 300s x 4, VIS 300s x 4, and NIR. (3s x 6) x 4. The spectra 
of bot h objects were reduced using the X-shooter pipeline 
vO.9.7 (IModigliani et alJfeoiQ ) following standard reduction 
procedures (flat-fields, wavelength calibration, etc.). The fi¬ 
nal wavelength calibration was carried out manually by iden¬ 
tifying several t elluric absorption line s using a solar spec¬ 
trum catalogue (IDelbouille et al.lll973l ). 

The UVB is corrected for instrumental response across 
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Figure 2. Ha line variability from three observation epochs for 
PDS 27: FORS2 (blue dashed line) and two X-Shooter epochs, X- 
shooter(a) (green dotted line) and X-shooter(b) (red solid line). 
For PDS 37, there is only one X-shooter data set (red solid line). 
A KC model (black dash-dotted line) of temperature 21000 K is 
also overplotted. The lines are shifted to the rest wavelength of 
the stars as in Fig. 1. 


the Balmer jump region (in order to measure the Balmer 
excess, which is a diagnostic for the strength of the accre¬ 
tion). To do this the observed spectra of the flux standards 
were divided through by the known flux from CALSPECQ- 
The resulting response curve was then applied to the tar¬ 
get of that night. To measure the Balmer excess, we do not 
need an absolute flux calibration, only the correct shape of 
the SED is needed. This method of calibration is therefore 
sufficient. We require only knowledge of the line strengths, 
profile shapes and calibrated wavelengths in the other arms. 
Therefore, this calibration is not needed in the other arms. 


3 RESULTS 
3.1 Spectral lines 

The spectra of both stars are remarkably similar showing 
numerous emission lines, while there is a paucity of pho- 
tospheric absorption lines. A representative set of lines is 
shown in Fig. |T] 

The strongest emission is found among the hydrogen 
recombination lines, Fe II, Ca II and [O I]. The Ha and H/3 
lines have P Cygni profiles, with absorption observed at ve¬ 
locities in excess of 500 kms -1 in the case of PDS 27, with 
some emission observed bluewards of the P Cygni absorp¬ 
tion. The Ha lines are very strong with line-to-continuum ra¬ 
tios of the order of 20 and an equivalent width (EW) of -70.4 
and -111.4 A for PDS 27 and PDS 37, respectively. Several 
strong Fe II emission lines can be seen in the visible portion 
of the spectra; the majority of these are double peaked in 
PDS 37, while they are all single peaked in PDS 27. Strong 
emission is also present in the Ca II triplet (8498, 8542, 


1 http://www.stsci.edu/hst/observatory/crds/calspec.html 


8662 A) and Ca II doublet (8912, 8927 A) in both stars. 
The Ca II triplet lines are blended with the Paschen lines at 
8502, 8545 and 8665 A. The Ca II line at 8542 and 8662 A in 
PDS 27 has a very weak blueshifted absorption. The other 
nearby Paschen lines at 8598, 8750 and 8863 A have simi¬ 
lar absorption components. Therefore, the absorption seen 
in the Ca II lines is more likely to belong to the Paschen 
series with which they are blended. He I lines at 5875, 6678, 
7065 and 10830 A are in absorption and blueshifted to the 
same velocities as observed in the hydrogen P Cygni pro¬ 
files. He I at 10830 A has a P Cygni profile in PDS 37 
with absorption observed up to -450 kms -1 , and the EW 
of the emission is slightly larger than that of the absorption. 
For PDS 27, He I at 10830 A has a very strong absorp¬ 
tion up to -700 kms -1 and is accompanied by a very weak 
redshifted emission component. As we go from low-order 
lines in the hydrogen Balmer series (Ha, H/3) towards the 
higher order hydrogen lines, we see a progressively weaker 
absorption component in each line. The only exception is 
that the absorption in Ha is less deep than in H /3, which we 
shall discuss later. The varying hydrogen line profiles pre¬ 
vent a reliable measurement of the stars’ radial velocity (14). 
Therefore, to obtain the systemic velocity, we measure the 
14 from the clean and symmetric Fe II and Ca II triplet lines 
in the X-shooter and FORS2 spectra. The measured 14 are 
corrected for heliocentric motion and local standard of rest 
motion for the dates of observation. These give local stan¬ 
dard of rest velocities ( Vlsr) of 47.2 ± 3.4, 45.3 ± 4.3, and 
63 ± 7 kms -1 for PDS 27 X-shooter(a), X-shooter(b) and 
FORS2, respectively. Both X-shooter(a) and X-shooter(b) 
values are in agreement with each other, while there is a dif¬ 
ference of ~ 18 ± 8 kms -1 between X-shooter and FORS2. 
For PDS 37, the Vlsr is found to be 11.7±5.5 kms -1 for X- 
shooter and 13.5±8.0 kms -1 for FORS2. For the final Vlsr 
to PDS 37 an average of the X-shooter and FORS2 values 
is taken and this results in a Vlsr of 12.6 ± 4.8 kms -1 . 

Both stars’ spectra are variable. An illustration of this 
is provided in Fig. [2]in which we show Ha at three observa¬ 
tion epochs for PDS 27 and two for PDS 37. The line profiles 
are variable with a deeper absorption component when the 
emission is weaker. For PDS 27 the Ha line-to-continuum ra¬ 
tio ranges from 13 to 18, and the EW increases from —70.4 
to —120.8 A from the X-shooter(b) to the FORS2 data sets 
(with the X-shooter(a) observations lying in between). For 
PDS 37 the line-to-continuum incre ases from 19.5 to 24 , and 
the EW from -111.4 to —122.6 A. IVieira et al.l (l201ll ) ob¬ 
served, in 1990, an EW of -88 and —105 A for PDS 27 and 
PDS 37 respectively. 

Variability of the He I line in PDS 27 also occurs be¬ 
tween all observations. Fig. [3] shows four of the He I lines, 
along with overplots from each observation epoch in order 
to highlight the variability. The He I lines 5875, 6678 and 
7065 A are present in both X-shooter and FORS2 spectra, 
while 10830 A is covered only by the X-shooter observations. 
The 5875 A line in the FORS2 data shows a broad blue 
absorption wing which is seen extending from -500 kms -1 
towards the red, but we cannot quantify how far it extends 
into the red due to the Na I lines blending with it. In the 
X-shooter(a) data set a clear absorption line is seen cen¬ 
tred at ~-150 kms -1 from the rest wavelength. The line in 
the X-shooter(b) is difficult to measure, due to the blending 
mentioned, but it appears to be centred at ~-100 kms -1 . 
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Table 1. The observed EW and the E(B-V) values obtained from 
the DIBs. Errors in the EW and E(B-V) are typically 5%. 


DIB 

(A) 

EW /E(B-V) 

(m A/mag.) 

PDS 27 

EW E(B- V) 
(mA) (mag) 

PDS 37 

EW E(B- V) 
(mA) (mag) 

5780.6 

579 

571.5 

0.99 

278.9 

0.48 

5797.1 

132 

190.7 

1.44 

221.6 

1.68 

5849.8 

48 

116.4 

2.42 

139.4 

2.90 

6089.8 

17 

24.1 

1.42 

23.4 

1.38 

6204.3 

189 

230.8 

1.22 

129.0 

0.68 

6013.7 

23 

242.0 

1.05 

206.5 

0.89 

6660.6 

51 

61.8 

1.21 

46.9 

0.90 

6993.2 

116 

131.7 

1.06 

100.0 

0.86 


3.2 Interstellar lines 






Figure 3. He I line variability is shown for PDS 27 of three ob¬ 
servation epochs: FORS2 (blue dashed line) and two X-shooter 
epochs, X-shooter(a) (green dotted line) and X-shooter(b) (red 
solid line). For PDS 37, there is only one X-shooter data set (red 
solid line) and one FORS2 observation (blue dashed line). A KC 
model (black dash-dotted line) of temperature 21000 K is over¬ 
plotted for the 5875, 6678 and 7065 A lines, while for the 10830 A 
line a black dash-dotted line shows just the continuum level. The 
lines are shifted to the rest wavelength of the stars as in Fig. 1. 

The 6678 and 7065 A lines show a similar behaviour in their 
velocities between the data sets, as demonstrated in Fig. 
[3] The 10830 A line only shows a small variability in the 
strength of the line. 

In the case of PDS 37, all the He I lines, present in both 
X-shooter and FORS2, show variability in the line centre of 
~50 kms -1 between them, with respective line centres being 
blueshifted from the central wavelength. Some of these show 
a very weak redshifted emission, which is seen to be much 
stronger in the 10830 A line. There is also strong variability 
in the line strengths, measured solely across the absorption 
component, where the strength is ~2.5 greater in FORS2 
than in X-shooter. 


The spectra of both targets contain interstellar absorption 
lines, most notably the K I doublet (7664.9, 7698.9 A), 
Na I doublet (5889.9, 5895.9 A) and the diffuse interstellar 
bands (DIBs); their origin comes from the material present 
in the interstellar medium (ISM). The strength of the DIBs 
increases with the amount of material, and therefore red - 
dening, along the line of sight Jjenniskens fe Desertl [l994l ). 
However, the D IBs can only provide a lower limit to the 
total reddening. lOudmaiier et all (1997) found for a sam¬ 
ple of young objects that the majority of the DIBs correlate 
with the reddening due to interstellar dust, but they re¬ 
main insensitive to the dust in circumstellar envelopes and 
their parental molecular clouds. Eight well-known DIBs are 
clearly identified in both objects; the measured EW for each 
line, along with the l ine streng th per unit E(B-V) deter¬ 
mined from Ijenniskens fe Desert] (Il994l ). and the resulting 
E(B-V) are given in Tabled] The overall E(B-V) is found 
to be 1.2 ± 0.16 for PDS 27 and 0.98 ± 0.37 for PDS 37. 
The spread of the values is fairly large with E(B-V) rang¬ 
ing, for PDS 27 and 37 respectively, from a low of 0.99 and 
0.48, for the 5780 A line, to a high of 2.42 and 2.9, for 
the 5849 A line. The very large colour excess derived from 
the 5849 A DIB is not en t irely unexpected as it was sug¬ 
gested by IChlewicki et al.l ( 1986 ) that the strength of this 
DIB could provide a good indication of the total reddening 
towards an object including the c ircumstellar component . 
This argument is strengthened by lOudmaiier et all (T997t ) 
who found in their study that this DIB was an exception¬ 
ally case which more accurately traced the total reddening 
for each star than any of the other DIBs. Therefore, in the 
overall E(B-V) reported above, the 5849 A DIB is not in¬ 
cluded. The Vlsr of the DIBs is measured to be 26.6 ± 5.3, 
27.3 ± 5.6 and 27.5 ± 2 kms -1 for PDS 27 X-shooter(a), X- 
shooter(b) and FORS2 spectra, respectively; all the values 
are in very good agreement with each other. For PDS 37, 
the Vlsr is found to be 0.5±2.8 and 6±4 kms -1 for the X- 
shooter and FORS2 spectra, respectively. This is consistent 
with the Vlsr observed in the sharp absorption lines from 
Na I 5889.9 and 5895.9 A, which is found to be 26.5 ± 1, 
26.8 ± 0.3 and 26 ± 2 kms -1 for PDS 27 X-shooter(a), X- 
shooter(b) and FORS2 spectra, respectively, and —4.5 ±0.5 
and —5 ± 1 for PDS 37 X-shooter and FORS2 spectra, re¬ 
spectively. 
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Table 2. The continuum polarisation of both objects measured in the following wavelength regions: B band centred at 4700 A, V band 
centred at 5500 A, R band centred at 7000 A, I band centred at 9000 A. The polarisation was measured over a wavelength range of 
250 A either side of the central wavelength, except for the B band where it was 125 A either side. All the errors are differential, where 
the systematic error is ~0.16% in polarisation and ~ 0.5° in angle (see the text in Section 2.1 for details). 


Object 

Pcont (%) 

B 

e° . 

cont 

Pcont(%) 

V 

e° . 

cont 

Pcont{%) 

R 

e° . 

cont 

Pcont (%) 

I 

e° . 

cont 

PDS 27 
PDS 37 

8.60 ± 0.03 
4.68 ± 0.04 

18.64 ± 0.11 
129.97 ± 0.27 

9.00 ± 0.01 
5.12 ± 0.01 

18.87 ± 0.01 
130.80 ± 0.01 

8.83 ± 0.01 
5.18 ± 0.01 

18.36 ± 0.04 
129.98 ± 0.06 

7.49 ± 0.01 
4.38 ± 0.01 

18.12 ± 0.03 
133.38 ± 0.05 


3.3 Continuum polarisation 


The continuum polarisation of the targets is shown in Fig. 
[|] In the figure, the spectropolarimetric data were rebinned 
using a coarse sampling, of ~ 250 A per bin, to minimize 
errors. It can be seen that the percentage of polarisation 
increases from short wavelengths, peaks in the V band and 
then decreases towards longer wavelengths. 

The observed continuum polarisation is a combination 
of interstellar polarisation and polarisation due to circum- 
stellar dust and electron scattering, if the dust and elec¬ 
trons deviate from circular symmetry. The nature of in- 
terstellar polari s ation is well described and explained by 
ISerkowski et alj (Il975l ). According to Serkowski’s law, the 
linear polarisation follows an empirical curve according to 
the following equation: 


PW 

p 

- 1 max 


= exp 



(3) 


where X ma x is the wavelength where the polarisation is at its 
maximum value, Pmax, and k is the width of the empirical 
curve. Typically, A max is ~ 5500 A, but can be in a range 
of 4500-8000 A. 

The Serkowski law provides an excellent fit to the data 
and is shown in Fig. [4] This combined with the large E(B- 
V) found towards these objects in Section 3.2 suggests that 
the bulk of the observed polarisation is indeed interstellar. 
The angle of polarisation changes slightly over the different 
bands; a change of 1 deg is observed for PDS 27 while a 5 
deg change is observed for PDS 37. The polarisation angle 
normally does not change as a function of wavelength for 
the ISM. The slight change in position angle with wavelength 
indicates that there is a contribution from another polarising 
agent to the total polarisation. The results of continuum 
polarisation are summarised in Tabled 


3.4 Line spectropolarimetry 

The polarisation spectra around the strongest lines, Ha, H/? 
and the calcium triplet are displayed in Fig. [5j There is 
a clear change in polarisation across the Ha line in both 
objects. The red peak and, where visible, the blue peak have 
almost the same polarisation value as the continuum, while 
there is a strong change across the absorptive component in 
both polarisation spectra and the ( Q,U) plots. The shift in 
the polarisation from continuum to the line is of the order 
of 1% for PDS 27 and 1.6% for PDS 37. 

Fig. [o] also shows that there is no clear change in the 
polarisation spectra across H/3 for both targets. Several dif¬ 
ferent ways of rebinning the data were tried, but no com¬ 



5000 6000 7000 8000 9000 

Wavelength (A) 


5000 6000 7000 8000 9000 

Wavelength (A) 


Figure 4. The polarisation and polarisation angles as a function 
of wavelength for PDS 27 and PDS 37. The spectropolarimetric 
data were coarsely sampled to improve the error bars. The black 
lines denote the best-fitting Serkowski law to the data. 


pelling evidence for line effects across H/3 could be found. 
Our explanation for this will be discussed in Section 5.1. 

The Ca II triplet has a lower ionisation potential than 
the Ha line, and as such it originates in a different region at 
a larger distance from the central stars. We did not detect 
any line effect across the Ca II triplet in PDS 27, whereas 
in PDS 37, the first two data sets show a possible line effect 
while the second two data sets do not show such a line effect 
(see Fig. [5}. Therefore, the current data are not sufficient to 
confirm whether there is a line effect or not across the Ca 
II triplet in PDS 37. We do not detect changes across any 
other lines in this polarimetric spectrum, which includes the 
Ca II doublet, He I, [O I] 6300 A and Fe II lines. 

The observed line characteristics of the targets in the 
FORS2 data are presented in Table [3] including an indica¬ 
tion whether a line effect is present. 
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Figure 5. The polarisation data of PDS 27 (upper half) and PDS 37 (lower half) are shown here. The data are presented as a combination 
of triplots (top) and ( Q,U) diagrams (bottom). In the triplot polarisation spectra, the Stokes intensity (I) is shown in the lowest panel, 
polarisation (%) in the centre, while the position angle (PA) is shown in the upper panel. The Q and U Stokes parameters are plotted 
against each other below each triplot. The data are rebinned to a constant value, as indicated at the top of each plot. The Ca II triplet 
in PDS 37 is presented in two colours: the first two data sets with possible line effect are shown in black solid line; the second two data 
sets with no line effect are shown in blue solid line. 
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Table 3. The line spectropolarimetry results: columns (2), (3), (5) and (6) list the Stokes (7) characteristics; 
columns (4) and (7) list line spectropolarimetry characteristics of each target. Errors in the EW measurements 
are typically 5%. 


Lines 

EW (A) 

PDS 27 
Line/cont. 

Line effect 

EW (A) 

PDS 37 
Line/cont. 

Line effect 

Ha 

-120.80 

18.0 

Yes 

-122.60 

19.4 

Yes 

H/3 

-14.56 

5.8 

No 

-12.27 

4.5 

No 

Ca 11(8498) 

-24.05 

6.0 

No 

-22.66 

4.4 

- 

Ca 11(8542) 

-23.97 

5.3 

No 

-25.43 

5.1 

- 

Ca 11(8662) 

-22.78 

5.2 

No 

-23.12 

4.8 

- 

Ca 11(8912) 

-1.11 

1.2 

No 

-0.81 

1.1 

No 

Ca 11(8927) 

-1.79 

1.3 

No 

-1.69 

1.3 

No 

He 1(6678) 

0.90 

0.9 

No 

1.40 

0.8 

No 

He 1(7065) 

0.27 

0.9 

No 

0.52 

0.9 

No 

O 1(6300) 

-0.36 

1.1 

No 

-0.46 

1.2 

No 

Fe 11(6456) 

-3.14 

1.8 

No 

-2.07 

1.5 

No 


4 ANALYSIS 

4.1 Stellar parameters 

4-1.1 Stellar velocity 

Due to the lack of photospheric absorption lines in the spec¬ 
tra of both stars, we used the clean and symmetric Fe II and 
Ca II triplet emission lines to estimate the systemic velocity. 

For PDS 27 we noted that in Section 3.1 both X-shooter 
observations provide a similar Vlsr, while there is a differ¬ 
ence of ~ 18 ± 8 kms -1 in Vlsr between X-shooter and 
F0RS2. To examine the reliability of this variation, it is im¬ 
portant to check the Ulsr of the interstellar lines as this 
should not vary over different epochs. In Section 3.2, we 
found a similar Vlsr of the DIBs in all epochs. These veloc¬ 
ities are consistent with the Vlsr observed in the sharp ab¬ 
sorption lines from Na I 5889.9 and 5895.9 A. The observed 
variation in Flsr of PDS 27 could be due to the star itself 
through changes in the circumstellar environment. On the 
other hand, it could be due to the presence of a second star 
in the system. Distinguishing conclusively between the two 
scenarios is challenging with only a few epochs of data, and 
requires additional observations. However, we found that the 
line profiles are almost symmetric in all epochs; hence, it is 
more likely that the changes seen are due to binary motions 
and are not due to the movements of circumstellar material 
around the star. Based on the above argument, we decide 
not to use our stellar velocities for PDS 27 for kinematic 
distance determinations. 

For PDS 37, the Vlsr is consistent within the errorbars 
for the X-shooter and FORS2 data. 


4-1.2 Distance 

Accurate stellar parameters are required for any quantita¬ 
tive analysis to be made about either age, accretion or evo¬ 
lution for these two enigmatic objects. First we will discuss 
the distance, which will later help constrain the total lumi¬ 
nosity. 

Using a Galactic rotation curve an estimate of the dis¬ 
tance to the stars can be made using the obser ved veloci- 
ties. In this work, we use the rotation curve from I Reid et al.l 
(120091 ). The observed Vlsr in PDS 27 is unlikely to trace the 


actual stellar velocity as argued in the previous section. For¬ 
tunately, molecular observations of 13 CO(1-0) and NH3(1,1) 
towa rds these objects are available dUrauhart et al.l 120071 . 
1201 lT) . The molecular lines in this case trace the parental 
clouds in a larger volume around the star compared to the 
optical lines. These authors observe velocities of 43.5 and 
5.4 kms -1 for PDS 27 and PDS 37, respectively. Using the 
molecular velocities the resulting kinematic distance for PDS 
27 is 3.17 (+0.66, -0.62) kpc, while for PDS 37 the distance is 
found to be 3.67 ± 0.95 kpc. For PDS 37 the stellar lines can 
also be used; these provide distances of 4.53 ± 0.77 and 4.74 
(+0.75, -0.90) kpc for the lines measured from X-shooter 
and FORS2, respectively. The molecular and spectral line 
measurements are in agreement with each other in PDS 37. 
For the final distance to PDS 37, an average of the three 
values is taken and this results in a distance of 4.31 ± 0.67 
kpc. 

An interesting consistency check is to see what the kine¬ 
matic distances to the DIBs along the line of sight towards 
the objects are. This is because the DIBs trace the mate¬ 
rial between us and the star, and should therefore be closer 
than the stellar and molecular distances. The observed Vlsr 
velocities of the DIBs in Section 3.2 provide kinematic dis¬ 
tances for PDS 27 of 1.76 (+0.58, -0.55) and 1.80 (+0.58, 
-0.55) kpc for the X-shooter and FORS2 data, respectively. 
For PDS 37 the kinematic distances of DIBs are found to be 
2.59 ± 1.4 and 3.76 ± 0.92 kpc for the X-shooter and FORS2 
data, respectively. These distances are indeed lower than the 
ones traced by the stellar and molecular lines, which is as 
expected. 

4-1.3 Temperature 

Spectral types of both stars are difficult to determine due to 
a general lack of photospheric absorption lines in the spec- 
tra. So far, the onl y spectral type given to these stars is by 
IVieira et al.l d2003lV . who assign them both a tentative B2 
classification based on their photometry. 

A rough lower limit on the t emperature can be found 
using the observed colour indices. IVieira et al.l d2003l i mea¬ 
sure a (B — V) 0 bs of 1.32 for PDS 27 and 1.52 for PDS 37. 
The DIBs provide colour excesses, E(B — V), of 1.20 ± 0.16 
and 0.98 ± 0.37 for PDS 27 and PDS 37, respectively. Com- 
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bining the observed and excess colours allows an intrinsic 
(B — V ) to be derived. This is only a lower limit as the DIBs 
trace the ISM, but not necessarily the total extinction. The 
highest ( B — V)i n t from the DIBs is found to be -0.04 in 
PDS 27 and 0.15 in PDS 37. This gives a lower limit to the 
temperatures which corresponds to B9-type star with a tem¬ 
perature of 10000 K and A5-type star with a temperature of 
8000 K for PDS 27 and PDS 37, respectively. The stars are 
likely to be even hotter than this, as this is a lower limit. 

An upper limit to the temperature can be placed using 
the Balmer jump in the spectra. This is because the Balmer 
jump decreases in size with increasing temperature (for tem¬ 
peratures greater than ~ 9000 K, where the Balmer jump 
is at its maximum). The observed size of the jump there¬ 
fore places an upper limit on the photospheric temperature. 
Moreover, a small jump can also be achieved by cooler stars 
combined with a flux excess due to accretion. The Balmer 
jump constraint provides an upper limit of 21000 K. A dis¬ 
cussion of the jump and possibility of an excess flux will be 
given in Section 4.2.1. 

The DIBs and Balmer jump constraints give an effective 
temperature range of 8000 K/10000 - 21000 K. In addition, 
the observed lines from the two objects suggest that they lie 
more towards the upper temperature limit. 

We detect N I emission lines at 8629 and 8683 A in both 
objects. These lines have a high excitation energy (~10 eV), 
and have only been obse rved in B-type HAeBe st ars, pre¬ 
dominately early types dHamann fe Perssonl I1992T) . Given 
that the latest spectral type in which these authors detected 
both NI emission lines was B 6 (corresponding to 14000 K), 
we infer a revised lower limit to the temperatures of PDS 27 
and PDS 37 of around 14000 K. 

Based on the above three arguments, which combine 
knowledge of photometry and spectra, we find the most 
plausible range to the temperatures to be 14000 - 21000 
K. We choose to adopt a temperature of 17500 ± 3500 K for 
both objects. 

In addition, the only absorption lines that are present 
in both spectra are from He I. If their origin were assumed 
to be from the photosphere, then their strength would peak 
at around 21000 K. However, the observed strengths are 
greater than predicted from stellar model atmospheres, and 
they are asymmetric. The blueshifted absorption extends all 
the way to -500 kms -1 and covers the same velocity range 
as the hydrogen lines. Therefore, an additional contribution 
to the absorption could be due to a strong wind. Although 
strong absorption has been regularly observed for T Tauri 
stars in the He I 10830 A line, such a strong and broad ab¬ 
sorption in a wind is not expected in the higher energy tran- 
sitio ns of He I in classical T Tauri stars (CTT s) by either the¬ 
ory (iKurosawa et al.ll20lH ') or observations dBeristain et al.l 
l200ll ') . This suggests that the stars are much hotter than 
CTTs, which is in agreement with the other temperature 
constraints adopted. 

4-1-4 Luminosity 

With a limit on both distance and temperature established 
for the stars, other stellar parameters can now be found. 
A start will be made with determining the luminosity us¬ 
ing the previously observed optical photometry. A total flux 
needs to be obtained, and due to the limited amount of pho¬ 


tometric observations an extrapolation into the UV must 
be performed to do this. To assess the UV contribution, 
Kurucz-Castelli (KC) models are used which provide a sur¬ 
face flux density. The B VRI photometry is then dereddened 
with this law until their shape matches the slope of the KC 
model for the adopted temperature. The U band is excluded 
due to possible effects of an excess due to accretion, while 
the longer wavelength photometry is also excluded due to 
being affected by an IR excess from the dust. To fit the 
models to the photometry, a scaling factor must be applied, 
as the models used are in surface flux while the photometry 
is observed flux. This scaling factor is just a ratio between 
the distance and radius squared, (D/R) 2 . For a temperature 
of 17500 ± 3500 K, a respective Av and D/R are found to 
be 5.03 ± 0.13 and 186.5 ± 21.0 for PDS 27, and 5.81 ± 
0.13 and 167.0 ± 18.9 for PDS 37. These reddening values 
are consistent with the DIB derived values, which were lower 
limits to the total reddening. Fig.^shows the result of dered¬ 
dening the photometry and a 17500 K KC model fit to the 
data. The luminosities are then calculated as a sum of the 
scaled KC model flux multiplied through by the distances 
found previously. Taking into account the errors on the dis¬ 
tance and temperature, the luminosities are calculated to be 
log(L*/LQ)= 4.39 ± 0.40 for PDS 27 and log(L«/L 0 )= 4.75 
± 0.39 for PDS 37. 

4-1-5 Location in the HR diagram 

With the luminosity established for each star, along with a 
well-defined temperature limit, the calculation of remaining 
stellar parameters through the use of PMS tracks can be 
performed. I n this work , two s ets of PMS tracks are adopted: 
one set by ISiess et~ak ( 2000 h. as they cover up to 12 M@, 
and the other set from lBernasconi fc Maedei] (|l996h . as they 
cover a larger luminosity and mass range (these are the 15 
and 25 Mq tracks in Fig.0. Along each mass track, there is 
a unique relation between the luminosity and temperature. 
Fig. [7] shows where PDS 27 and PDS 37 lie on the PMS 
tracks. Based on their position, we find a mass of 15.3 (+5.4, 
-4.4) M@ for PDS 27 and a mass of 21.1 (+11.0, -5.3) M@ 
for PDS 37. The radius is found straightforwardly from the 
D/R values, determined in the previous section, by using the 
distance. This provides 17.0 ± 4.0 R@ for PDS 27 and 25.8 ± 
5.0 R@ for PDS 37. With a mass and radius determined, the 
log(g) is calculated to be 3.16 ± 0.27 for PDS 27 and 2.94 
± 0.35 for PDS 37. The errors in log(g) are dominated by 
the errors in distance. The position of the stars in the HR 
diagram suggests that they are very young objects, which 
may be evolving into massive O-type stars. 


4.2 Accretion rates 

Measurement of the accretion rates is done by two separate 
methods: (1) a direct relationship to the measured Balmer 
excess and ( 2 ) using the strength of accretion-tracing lines. 

4-2.1 Accretion rates 1: Balmer excess 

The first method of obtaining accretion rates is to mea¬ 
sure the Balmer excess, A Db ■ This is defined as the ex- 
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Figure 6. An example of the SED integration for extracting a dereddened total flux for the case of PDS 27 (left) and PDS 37 (right). 
The red diamonds are the observed photometry, the blue crosses are the dereddened photometry and the solid black line is a KC model 
of 17500 K. 


Table 4. Lines tracing accretion. Errors in the EW measurements are typically 5%. Errors in the line flux are ~80%, for both stars, due 
to the large temperature and distance errors assigned. EW„j, s , EW, rt i and EW cor are observed, intrinsic photospheric absorption and 
corrected EW, respectively. 


Line 

Star 

Spectra 

EW 0 ), S 

(A) 

EW, nI 

(A) 

EWcor 

(A) 

E line 

(W/m 2 ) 

log(L line) 

(L©) 

log (Z/acc) 

(L©) 

log(Macc) 

(M©/yr) 

Ha 

PDS 27 

X-shooter(a) 

-75.20 

4.20 

-79.40 

l.lxlO -13 

1.54 ±0.10 

3.96±0.07 

-3.49±0.20 

Ha 

PDS 27 

X-shooter(b) 

-70.40 

4.20 

-74.60 

l.OxlO -13 

1.52±0.10 

3.93±0.07 

-3.52±0.20 

Ha 

PDS 27 

FORS2 

-120.80 

4.20 

-125.24 

1 . 8 xl 0 -13 

1.74±0.10 

4.18±0.06 

-3.27±0.19 

O 

05 

o 

o 

PDS 27 

X-shooter (a) 

-0.28 

0.00 

-0.28 

4.6xl0 -16 

-0.84±0.10 

3.85±0.13 

-3.60±0.23 

[OIJ6300 

PDS 27 

X-shooter(b) 

-0.29 

0.00 

-0.29 

4.8xl0 -16 

-0.82±0.10 

3.87±0.14 

-3.58±0.23 

[OIJ6300 

PDS 27 

FORS2 

-0.36 

0.00 

-0.36 

5.9xl0 “ 16 

-0.73±0.10 

3.97±0.16 

-3.48±0.24 

Bry 

PDS 27 

X-shooter(a) 

-3.80 

4.72 

-4.22 

8.3xl0 -16 

-0.58±0.10 

3.02±0.16 

-4.43±0.24 

Bry 

PDS 27 

X-shooter(b) 

-3.80 

4.72 

-4.22 

8.3xl0 -16 

-0.58±0.10 

3.02±0.16 

-4.43±0.24 

Ha 

PDS 37 

X-shooter 

-111.40 

4.13 

-115.53 

2 . 0 xl 0 -13 

2.07T0.08 

4.54±0.04 

-2.87±0.22 

Ha 

PDS 37 

FORS2 

-122.60 

4.13 

-126.73 

2 . 2 xl 0 -13 

2.11 ±0.08 

4.58±0.04 

-2.83±0.22 

[OI]6300 

PDS 37 

X-shooter 

-0.47 

0.00 

-0.47 

9.6xl0 -16 

-0.25±0.08 

4.51±0.36 

-2.89±0.42 

[OIJ6300 

PDS 37 

FORS2 

-0.46 

0.00 

-0.46 

9.4xl0 -16 

-0.26±0.08 

4.50±0.35 

-2.91±0.41 

Br 7 

PDS 37 

X-shooter 

-8.50 

4.64 

-9.28 

1 . 2 xl 0 -15 

-0.15±0.08 

3.41±0.49 

-4.00±0.53 


cess of energy seen across the Balmer ju mp region, and 
it ha s been observed in numerous CTTs (IGulIbring et al.l 


2000l: Herczee & Hillenbrand 20081: Inalebv et al. 

2013) 

and HAeBe stars ( 

Garrison 19781: iMuzerolle et al.l 

2004; 

Donehew & Brittain 

2011: Mendieutfa et al. 2011. 120141). 


The excess energy is often explained to be caused by mag- 
netospheric accretion (MA), where disc material is fun¬ 
nelled via accretion columns on to the star where they 
shock the photosphere, with the final result of this being 
a hotspot of increased te mperature that peaks in the UV 
dCalvet fe Gullbrinelll998l ). Therefore, A Db serves as a di¬ 
rect tracer of the accretion rate. 

The method of measurement here follows 
iDonehew fc Brittainl d201ll ). where the observed spec¬ 
tra and a model atmosphere of the star’s temperature 
are normalised to 4000 A, and then the spectra are 
adjusted by a reddening law until the slope of the SED 


matches that of the model between 4000 and 4600 A; this 
accounts for any reddening in the s pectra, and has been 
shown to be reddening ind ependent dMuzerolIe et al.ll2004l : 
IDonehew fc Brittainl l201ll ). Finally, the Balmer excess is 
measured by the equation: 

A D b = 2.51og (4) 

where F * is the flux of the photosphere at 3600 A; F, a CC is the 
flux produced due to accretion, at 3600 A, such that F„ ±F acc 
is the observed flux. Fig. [8] shows the Balmer jump region 
of each object and the expected intrinsic photosphere, given 
by a KC model. The spectra is binned to 5 A resolution in 
order to clearly see and measure the excess, and the region 
measured is from 3605 to 3665 A to avoid the poor SNR 
region where the echclle orders crossover. A Db is measured 
to be 0.17 ± 0.06 for PDS 27 and 0.16 ± 0.04 for PDS 
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Figure 8 . The Balmer jump region for PDS 27 (left) and PDS 37 (right). The targets are shown as the green solid line while the magenta 
dotted line is a KC model of 17500 K. The plot also highlights the extreme emission line nature of the two objects. The error in the 
Balmer excess displayed is solely the error due to the measurement between the two spectra. 
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Figure 7. Show s the P MS tracks of lBernasconi fc Maedej (119961) 
and ISiess et alJ (|2000l ). along with the positions of PDS 27 in 
magenta dashed line and PDS 37 in green solid line. The dotted 
lines show where the objects lie for the determined temperature 
range, while the error bars show the allowed changes in luminosity 
given the errors in distance. 


37. However, because of the temperature range adopted, the 
measured A Db can vary by more than this error (see Fig. n. 
Therefore, we assign larger, more generous, errors to these 
measurements of 0.17 (+0.13, -0.16) for PDS 27 and 0.16 
(+0.22, -0.15) for PDS 37. 

As mentioned in Section 4.1.3, an upper limit was 
placed on the temperature of the stars by considering the 
observed size of the Balmer jump region. For the case of 
stars with an excess flux in this region, due to accretion, we 
can instead measure the Balmer excess, which is the amount 



14000 16000 18000 20000 22000 24000 

Temperature (K) 


Figure 9. The measured Balmer excess changes for each star 
depending upon the adopted temperature. PDS 27 is shown in 
magenta dashed line, and PDS 37 in green solid line. For both 
stars, a zero in Balmer excess is achieved around 21000 K. 


the Balmer jump is filled in by the excess flux. In Fig. 0the 
observed Balmer excess for each star is shown for different 
temperatures (i.e. different Balmer jump sizes) to highlight 










































































































































12 K. M. Ababakr, J. R. Fairlamb, R. D. Oudmaijer and M.E. van den Ancker 


that there is a limit in temperature beyond which no excess 
can be observed. 

The relationship between A Db and M ac c depends upon 
the stell ar parameters of the sta r, and a full description is 
given in iMendigutia et ahl (l201ll l . The basics are that the 
derivation of M acc from A Db is unique for each combina¬ 
tion of stellar parameters, i.e. in a 20000 K star the same 
measured A Db will yield a much higher accretion rate than 
a 10000 K star because it requires a much higher accretion 
flux for the same visible A Db- Because of this, it also means 
that for hotter objects smaller A Db values are expected due 
to the sheer amount of accretion flux required to be notice¬ 
able above the photospheric level. 

For the parameters determined in this work, the 
log(M acc ) found are -3.96 (+0.76, -1.32) for PDS 27 and 
-3.56 (+0.60, -1.62) for PDS 37. Filling factors are also de¬ 
termined of f = 0.40 (+0.55, -0.40) and 0.40 (0.54, -0.39), 
for PDS 27 and 37, respectively, when using a typical in- 
ward flux of accretion column energy of 10 12 erg/cm 2 /s 
(IMendigutia et aj]|201 il l . These are large but acceptable fill¬ 
ing factors. 


4-2.2 Accretion rates 2: line luminosity 

For the second method, there are established relationships 
between accretion rate and emission line strength. This 
is because during the accretion process large amounts of 
energy are released which can excite the circumstellar 
environment. This excitation a nd excess of emission in 
CTTs can be explained by MA (ICalvet fe Gullbrindll998l : 
iMuzerolle. Calvet fe Hartmannl 1 1 9981 ). Well-established 


line stre ngth relationships have been found for brow n 


dwarfs et al. 20041. 200f& Wh 

ite & 

3asril 

200311. 

CTTs dMuzerolle. Hartmann & Calvet] 

19§1; 

Dahm 

20081: 

llnglebv et al.l 20131: lAlcala et al.l 201411 

and HAeBe stars 

(Garcia LoDez et alJ 120061: IMendigutia 

et al.! 1201 ill. 

Since 


the two objects in this work are ear l y-type Herbig stars, 
the relationships of IMendigutia et ahl (l201ll l are adopted. 
However, it should be noted that these relationships may 
not hold for the targets in this work which are highly 
luminous compared to the stars from which the relationship 
was derived. This is an aspect we aim to test with the direct 
A Db derived accretion rate. 

To obtain the line luminosities, a dereddened line flux is 
needed. This requires a line strength to be measured, which 
must be corrected for intrinsic line strength. The observed, 
intrinsic and corrected line strengths are given in Table [4] 
For the Bry line, it should be noted that the intrinsic cor¬ 
rection must also take into account the IR. excess in the K 
band. This is done following the steps in lGarcia Lopez et al.l 
(l2006li . The flux density is then calculated from an inter¬ 
polation between the two closest photometric points; with 
the exception of the Br 7 line for which no interpolation is 
made, instead the value of the A-band magnitude is used. 
The photometry is dereddened by the Av values found in 
Section 4.1.4 prior to this step. The line flux is obtained 
by multiplication of the interpolated flux density and the 
line strength, and the line luminosity is then achieved by 
multiplying the line flux by 47 tD 2 . 

Once the line luminosity is established, the accretion 
luminosity for each line is extracted via the relationships of 
IMendigutia et ahl d201ll l which take the form of log(L acc ) = 


A + B log(Liine), where A and B are constants specific to 
each line. Finally, the accretion rate is calculated via the 
equation: M acc = T aC c-R*/GM*, where Gis the gravitational 
constant and M acc is the accretion rate. Table[4]gives details 
of the values used to derive the accretion rates. The accretion 
rates derived from the Ha, [O I] and Bry lines are in the 
range of M acc ~ 10 -3 10 -4 ' 5 Mg/yr. 


5 DISCUSSION 

We have presented a spectroscopic and spectropolarimet- 
ric study of two early-type Herbig stars, and found that the 
spectra of both objects are remarkably similar, strongly sug¬ 
gesting that both are in the same evolutionary phase. The 
spectra are dominated by strong emission lines, and do not 
reveal any photospheric absorption lines. Indeed, the only 
lines that are fully in absorption are the interstellar DIB 
lines. 

The picture that emerges is that both objects undergo 
strong winds as attested by the clear, and numerous, P 
Cygni profiles observed. P Cygni absorption extends in both 
stars to hundreds of kms -1 . The He I line at 10830 A is the 
only helium line seen as a P Cygni profile. The remaining He 
I lines at 5875, 6678 and 7065 A are present only in absorp¬ 
tion (see Fig. ®. These lines are broader than the expected 
intrinsic lines from the photosphere. If they originate from 
the wind, they must probe the very hot inner wind. 


5.1 Circumstellar geometry 


The majority of the observed Fe II lines are double peaked 
in PDS 37, while they are single peaked in PDS 27 (see Sec¬ 
tion 3.1). This suggests that PDS 37 has an inclined disc 
that we are viewing close to edge-on. This interpretation 
is supported b y a high inclination angle being reported by 
lllee et all (120131 '). based on first overtone CO emission. In 
addition, the continuum polarisation shows a change in posi¬ 
tion angle with wavelength. The observed polarisation must 
have an intrinsic component as normally the interstellar po¬ 
larisation angle is constant with wavelength. This change is 
stronger in PDS 37 than PDS 27 as illustrated in Fig. [I] Fur¬ 
ther information on the geometry of the circumstellar envi¬ 
ronment can be achieved by performing spectropolarimetry 
across the lines. 

Spectropolarimetry is very effective in revealing asym¬ 
metries in the scattering material. This technique is based 
on the concept that free electrons very close to the stars, 
where the density of free electrons is high due to hydrogen 
ionisation, scatter the continuum radiation from the stellar 
photosphere and polarise it. The electron scattering reg ion 
is of the order of a few stellar radii (ICassinelli et al.lll987f ). In 
the case of a face-on disc, or if surrounding material is spher¬ 
ically symmetric, no net polarisation can be detected, while 
in the case of inclined flattened structures a net polarisation 
can be observed. The emission lines are less polarised than 
the continuum radiation because they face fewer electrons 
than the continuum and the ionised region is larger than 
the electron scattering region. Clarke fe McLeanl rtl974T ) and 
IPoeckert fe Marlborough! (Il976h first used spectropolarime- 
try in a study of Be stars, which allowed them to infer that 
these stars are surrounded by discs. 
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A change in polarisation across the absorptive com¬ 
ponent of the emission line can not be explained by the 
above concept, and is often suggested to be the result of the 
McLean effect. In brief, the absorption takes the direct, un¬ 
scattered, light from the star away from the beam; however, 
the flux does not necessarily reach zero. Some of the light 
emitted in other directions can be scattered into our line of 
sight by free electrons, provided the electrons are in a non- 
spherical distribution. Therefore, the observed light in the 
absorption reg ion will be more polarised than the continuum 
jMcLeanlll979r ) . There are many objects reported a s having 
a McLean effect; see for example I Vink et al.l (l2005h 

Both PDS 27 and PDS 37 show a strong McLean-type 
effect across the absorptive component of Ha (see Fig. 0 
and the change is as broad as the absorption. This indicates 
the presence of intrinsic polarisation which in turn implies 
that the ionised material around the star is asymmetrically 
distributed. The strength of the line effect in Ha is of the 
order of 1% in PDS 27 and 1.6% in PDS 37. This supports 
the ideas that PDS 37 has a more inclined disc than PDS 27, 
as we expect more intrinsic polarisation in a more inclined 
disc. The slope of the loop from the continuum to the line in 
the ( Q,U) diagram represents the intrinsic polarisation for 
the continuum polarisation (see Fig. 0. The intrinsic con¬ 
tinuum polarisation angle represents the electron scattering 
angle in the circumstellar medium. We determine an intrin¬ 
sic polarisation angle of ~ 77° for PDS 27 and ~ 56° for 
PDS 37. The angle has been measured from the continuum 
to the line in the ( Q,U) diag ram in Fig. 0 [see the discussion 
in I Wheelwright et alj (120111 )]. 

We would expect a similar line effect in H/3 as in Ha, 
since it arises from approximately the same volume. How¬ 
ever, we do not see any changes across the line profile, nei¬ 
ther in absorption nor in the emission (see Fig. 0 . A point 
noteworthy is that the absorption component in H/3 is deeper 
than that of Ha, which is not intrinsically expected as the 
latter has a greater optical depth than the former. 

In the present cases, we are able to explain this rela¬ 
tively straightforwardly, by using the information that the 
strong emission from Ha is scattered into our line of sight, 
filling in its P Cygni absorption component with scattered, 
polarised emission. In contrast, H f5 is a much weaker emis¬ 
sion line, and therefore far fewer photons are available to be 
scattered into the line of sight compared to Ha. As a con¬ 
sequence, weaker, if any, changes in polarisation across the 
P Cygni absorption would be predicted in this model. This 
explains a non-detection of the line effect across the absorp¬ 
tion component in H/3. This argument can also explain why 
the Ha line effect in PDS 37 is more obvious than the one 
in PDS 27, as the latter line has deeper absorption com¬ 
ponents. For all of the other emission lines covered by the 
FORS2 spectropolarimetry of both stars, we do not detect 
a line effect. 

To summarise, we do not see a line effect across the H/3 
line as the majority of the photons in the absorption compo¬ 
nent are unscattered light. The observed He I emission lines 
are weak and they do not show any line effect, and the line 
effect across the Ca II triplet in PDS 37 is not compelling. 
However, from the characteristics of the observed spectral 
lines, along with the continuum polarisation and the line 
spectropolarimetry across the Ha line, we can confirm that 
these two objects have an asymmetric circumstellar environ¬ 


ment. We identify this to be a circumstellar disc, whose size 
is at least of the order of several stellar radii. 

5.2 Stellar properties 

Based on their position in the HR diagram, we can infer 
that both of the targets are young objects with an active 
circumstellar environment. This is evidenced by high accre¬ 
tion rates of (Macc) = 1CV 3 -10 -4 ' 5 Mgyr -1 ; strong winds, 
with P Cygni absorption of around 90 and 80% of continuum 
level for PDS 27 and PDS 37, respectively, in H/3, and also 
maximum respective velocities of ~-700 and ^-400 kms -1 
measured in the Ha line; and variability in line strengths 
and velocities seen in both the hydrogen and the He I lines. 
In particular, PDS 27 has shown unusual variability in the 
He I lines where the velocity of the line centre and the line 
strength change dramatically between epochs. This type of 
variability in line strength has be e n obs erved previously in 
other objects by iKucerova et ahl (l2013l h but an origin of 
the varying line strength cannot be reached with only three 
epochs of data. PDS 27 also shows variability in the observed 
Flsr between the X-shooter and FORS2 spectra. This could 
be due to the presence of a secondary companion. However, 
with only three epochs of data, we cannot establish con¬ 
clusively the binary nature of PDS 27; additional data are 
required to determine the frequency and magnitude of the 
velocity variations. Unlike PDS 27, PDS 37 did not show 
any variability in the observed Ulsr between X-shooter and 
FORS2 spectra. 

Their young age is further supported by their large radii 
and low surface gravities. These properties are typical of 
very young objects where their puffed-up envelope is still 
undergoing contraction a nd will shrink as i t progresses to¬ 
wards the main-sequence dDavies et al.ll201(il j. The accretion 
rates obtained from the Ha and [O I] 6300 A are found to 
be the same across all epochs, within the errors, for both 
PDS 27 and PDS 37. These two lines also show that PDS 
37 has an accretion rate ~0.6 dex greater than PDS 27, 
which is comparable to the difference seen between the ac¬ 
cretion rates derived from the Balmer excess of ~0.4; both 
are within the errors. The Bry line does not follow the same 
trends as it shows a decrease in accretion rate, in both ob¬ 
jects, over both the other two lines and the Balmer excess 
measurement (see Table 0 for all of the emission line derived 
accretion rates). 

The accretion rates derived using the Bry line are 
higher than any other reported accretion rates in HBe stars, 
which also use the Bry line, by ~l-3 orders of mag nitude 
jGarcia Lopez et al.l 120061 : iDonehew fc Brittainl 1201 lh . The 
Balmer excess method results in log(M acc ) =-3.96 (+0.76, 
-1.32) for PDS 27 and -3.56 (+0.60, -1.62) for PDS 37. These 
values agree with the respective accretion rates derived from 
all of the lines analysed. The Bry line in both objects shows 
a difference of an order of magnftude less in accretion rate 
between the other lines. A possible breakdown in the rela¬ 
tionship between line l uminosity and accret i on lum inosity 
has been suggested by IDonehew fc Brittainl (1201 il l , where 
the luminosity increases towards the early-type HAeBe star 
regime. However, the two stars presented in this work still 
follow this relationship. We cannot conclude with just these 
two objects whether the relationship holds or not, and such 
is beyond the scope of this paper. Instead, PDS 27 and PDS 
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37 will form part of an upcoming paper exploring these rela¬ 
tionships in a larger sample (Fairlamb et al., in preparation). 


6 CONCLUSIONS 

This work presents new data on two previously poorly stud¬ 
ied HAeBe candidate stars. By combining X-shooter spec¬ 
troscopy and FORS2 spectropolarimetry, our understanding 
of properties in the circumstellar environments and the stel¬ 
lar parameters of the stars themselves has improved. The 
main findings in this work are as follows: 

• The targets are both very hot 17500 ± 3500 K with 
a spectral type of around B2. When combined with their 
high luminosities, of log(L,/LQ) ~ 4.5, they are placed in 
a secluded area on the HR diagram where very young ob¬ 
jects lie. They are swollen objects with large radii, 17.0 ± 
4.0 and 25.8 ± 5.0 Rq, and large masses, 15.3 (+5.4, -4.4) 
and 21.1 (+11.0, -5.3) Mq, for PDS 27 and 37 respectively. 
Their young nature is supported by the high accretion rates 
calculated here, which are among the highest measured in 
HAeBe stars to date, M a cc ~ 10 -3 10 -4 ' 5 MgyU 1 . 

• Spectropolarimetric line effects are detected in the Ha 
line for both objects. This polarisation property is explained 
by the McLean effect; line photons are scattered into the 
line of sight, increasing the polarisation across the P Cygni 
absorption. The spectropolarimetry indicates that the cir¬ 
cumstellar environment is of a flattened structure, likely a 
disc. 

• No line effects are seen across lines other than Ha. This 
is consistent with the notion that the spectropolarimetric 
behaviour in both cases is due to the McLean effect. The 
effect requires a strong emission line filling in absorption 
and a line with great optical depth so that the scattered 
photons will dominate the polarisation. It appears that in 
the cases studied here, only Ha fulfils both requirements. 

• Line variability is seen in both objects, which we as¬ 
cribe to an active circumstellar environment. Extreme line 
variability of the He I species is seen in PDS 27. In addition, 
PDS 27 also shows a variability in the observed Vlsr be¬ 
tween X-shooter and FORS2; this could be due to the fact 
that the system is binary. 

Finally, we have presented two young PMS stars, whose 
spectral classification as early B-type is consistent with them 
belonging to the group of HAeBe stars. Their masses and lo¬ 
cations in the HR diagram indicate that they will continue to 
evolve to become O-type stars. As they are optically bright, 
these stars provide important test beds to study the forma¬ 
tion of the most massive stars. 
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